Character of Atomic Vibrations in a Lennard-Jones Glass 
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Lennard-Jones glasses (made on a computer by 
quenching from liquid state coordinates) are studied 
in harmonic approximation. Vibrational eigenfrequen- 
cies and eigenvectors are found by exact diagonaliza- 
tion for models with periodic boundaries and N=(500, 
2048, and 6980) atoms. We analyze the density of 
states, mobility edge, and bond-stretching character 
of the normal modes. In agreement with older work of 
Grest, Nagel, and Rahman, the upper 7% of the modes 
are localized, and the rest, delocalized. The modes can 
not be differentiated by any property or quantum num- 
ber except their eigenfrequency. More specifically, in 
a given narrow frequency interval, all modes are glob- 
ally identical. Transverse or longitudinal character, for 
example, disappears. 
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I. INTRODUCTION 

Vibrational properties of disordered media have 
been reviewed by various authors [1-4]. Subse- 
quent to most of these reviews, Grest, Nagel, and 
Rahman (GNR) [5-10] did a careful characteriza- 
tion of the normal modes of the Lennard-Jones 
(LJ) glass. Our previous numerical studies, mo- 
tivated by heat conductivity experiments and by 
theoretical controversy, focussed on vibrations in 
amorphous Si (Am-Si) [11]. Here we present re- 
sults for a very large LJ glass model, with local 
atomic structure very different from Am-Si. These 
results confirm the GNR picture. We use the term 
"glass" simply to mean an amorphous solid, and 
the abbreviation am-LJ to denote the LJ glass. 

The model studied is a monatomic 6-12 
Lennard-Jones system, V(r) — 4e[(er/r) 12 — 
(cr/r) 6 ], with potential parameters suitable for Ar- 
gon: e/k B = 125.2 K and a = 0.3405 nm. This po- 
tential gives a satisfactory description of crystalline 
Ar [12]. In its face-centered cubic (fee) structure, 
each atom is surrounded by 12 nearest neighbors 
located a little inside the minimum of the V(r) at 



r = 2 1 / 6 er. The second shell is more distant by 
2 i/2 

Standard molecular dynamics simulations were 
done on systems with N=500, 2048 and 6980 
atoms, enclosed in a cubic box with periodic 
boundary conditions. A truncated (r c = 2.6a) 
LJ potential was used. The simulated density was 
p = 42 x 10~ 3 mol/cm 3 , corresponding to the re- 
duced value p*—pa 3 ~ 1. After an equilibration 
run at T = 400 K using NVE molecular dynamics, 
gradual cooling was performed down to a temper- 
ature (T = 150 K) slightly above the melting tem- 
perature. Finally, conjugate gradient minimization 
found the glassy minimum configuration. 

The resulting structure of am-LJ is somewhat 
like that of a smeared-out argon crystal. The den- 
sity is about 5% smaller. There may be icosahedral 
aspects in the local packing [13], but our models 
are more similar to the "LJ glass" of Simdyankin 
et al. [14] than to their "IC glass" with strongly 
icosahedral features. A shell of radius r contains 
on average (4ir/3)pr 3 atoms, where p is the num- 
ber density. The mean distance from an atom to 
its iV'th ncig hbor is thus (SN/Airp) 1 / 3 . 
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FIG. 1. mean distance to the A^h neighbor, for 
am-LJ (solid curve) and the corresponding crystal 



(dashed lines). Error bars give the variance < 5r > 
in the distribution of TV'th neighbor positions . 
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For am-LJ we calculated the mean distance to 
the JV'th neighbor, and the variance < 8r 2 > 1 / 2 of 
this quantity, which are plotted in Fig. (1). The 
closest 10 neighbors cluster in a narrow interval 
near the distance 3.7A found in crystalline Ar. The 
11th to 13th neighbors have less well-defined posi- 
tions. The 15th to 18th atoms have the least well- 
defined positions, lying near the second neighbor 
shell of the crystal. More distant neighbors start 
to be packed as in a random gas, with distance 
starting to scale like TV 1 / 3 . 

Normal modes of vibration were calculated in 
harmonic approximation, using the same periodic 
boundary conditions used to compute the posi- 
tions. The calculations were done on Njal, a Be- 
owulf cluster with 81 dual pentium III processors 
running at 0.933 GHz. The Scalapack routine 
PDSYEV was used. The largest system size, 6,980 
atoms, corresponds to a real-symmetric dynami- 
cal matrix of dimension 20,940, the largest that 
fits conveniently into Njal's memory. This matrix 
required about 100 minutes to diagonalize. The 
vibrational density of states is shown in Fig. (2). 
All three models give essentially the same result, 
which is shown here for the 6980-atom model, and 
was first found by Rahman et al. [15]. The re- 
sults for the glass are compared with the density 
of states of the crystal, which closely resembles ex- 
periment on crystalline Ar [16]. However, the LJ 
glass has a surprisingly different density of states, 
being featureless, with no distinction between lon- 
gitudinal and transverse modes. 
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FIG. 2. Density of vibrational eigenstates of am-LJ 
and the corresponding fee Lennard- Jones (argon-like) 
crystal. The dashed vertical line indicates the position 
of the mobility edge in am-LJ. 

To pursue this issue, it is interesting to look 
at the "bond-stretching parameter" S a of the a'th 
normal mode, defined as [17] 



a L £<^>i«-^)i 2 J 

where uf is the displacement eigenvector on the 
atom with position label i and position Ri. The 
unit vector hij points in the direction Ri — Rj . In 
other words, S a measures for mode a the relative 
amount of stretching of each near-neighbor bond. 
For the crystal, this is rigorously defined by using 
the 12 bonds from each atom to its first neighbors. 
For the glass, we sum all neighbors whose separa- 
tion is less than the mean distance to the 12.5th 
neighbor. The results are plotted in Fig. (3). 
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FIG. 3. Bond-stretching character versus frequency 
for the normal modes of the glass (narrow dark band) 
compared with the bond-stretching character of the 
normal modes of the crystal (broad light band). 

This figure shows a very striking difference be- 
tween the crystal and the glass. Crystalline normal 
modes can be labelled by the wavevector k. In Ar, 
there are only three normal modes at each k. They 
have usually quite different bond-stretching char- 
acters, and can normally be classified as longitu- 
dinal (L) or transverse (T). This leads to a highly 
differentiated spectrum of states, seen in the figure 
as a broad range of values of S a at each frequency. 
In the glass, the bond-stretching character S a col- 
lapses to a narrow interval whose width diminishes 
to zero as N increases to infinity. Lacking the 
quantum number k, each mode has totally mixed 
character. This property was noticed previously 
for other glasses [18,19,11]. In an infinite glass, 
with infinitely many modes at each frequency, one 
could make a unitary rotation (within the subspace 
of fixed frequency) to a basis where certain modes 
would look approximately pure L in some local re- 
gion. However, at long enough distances this char- 
acter is completely lost. At low frequencies, puri- 
fied modes defined this way could be chosen have 
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a fairly sharply defined local value of k. At large 
enough distances, the mode would evolve into a 
spherical mixture of different directions of k. As 
frequency increases, there is a crossover (the "Ioffe- 
Regel" (IR) crossover) to the regime where the lo- 
cal k can hardly be defined because it decays in a 
distance of one or two neighbors. The IR crossover 
and other low-frequency properties will be the sub- 
ject of a future study. 

It was noticed by GNR that a given normal 
mode, insofar as it was L, expressed its L character 
at specific wavevectors, and insofar as it was T, ex- 
pressed that character at other wavevectors. This 
enabled them to draw pseudo-dispersion curves 
showing very different character for L and T as- 
pects. This does not contradict our view that a 
given mode cannot be classified. Each mode is si- 
multaneously L and T, probably even at the same 
point in space unless purified as described above. 
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FIG. 4. Inverse participation ratio on logarithmic 
scale versus eigenfrequency for all normal modes of sys- 
tems of various sizes. 

In principle, normal modes of vibration, just 
like electron eigenstates in quantum theory, can 
be rigorously classified as extended (E) or local- 
ized (L). In three dimensions the vibrational spec- 
trum has sharp E/L boundary frequencies ("mo- 
bility edges," u c ) separating these two kinds of 
modes. The IR crossover is low in the spectrum 
of a glass, but the E/L boundary is usually high in 
the spectrum. 

An hypothesis that the mobility edge coincides 
with the IR crossover [20] has not proved true for 
any model that has been studied numerically, but 
has remained a popular misconception. Our large 
system provides a useful confirmation that pre- 
vious studies have correctly located the mobility 
edge. We compute the participation ratio p a for 
each normal mode, defined as 



The values of p a are plotted on a logarithmic scale 
versus eigenfrequency for all three models in Fig. 
(4). 

Notice first that the value of p a hovers at slightly 
higher than 0.57V for a significant part of the spec- 
trum (frequencies 2.5 to 6.5 meV). These are very 
well delocalized states whose amplitude does not 
vary too much from atom to atom. At low fre- 
quencies (below 1 meV for the 6980-atom model, 
and pushing up to 2 meV for the 500-atom model) 
there is a trend to lower values of p a , and there oc- 
cur some idiosyncratic modes with p a reduced by 
10 or more. Their interpretation remains difficult 
[21]. At least some of this behavior has to be at- 
tributed to finite-size artifacts. We hope to make 
a detailed study later. Finally, at high frequen- 
cies, very low values of p a occur. These states are 
strongly Anderson localized. The mobility edge is 
estimated to be at uj c = 8.2 ± 0.1 meV. This es- 
timate was found by fitting a straight line to the 
logarithmic data at the point where p a w 0.257V. 
Extrapolating this line back to the horizontal line 
p a = 0.57V, the intersection occurs at around 8.2 
meV, with uncertainty which diminishes as the sys- 
tem size increases. This value of u c corresponds to 
having the lower 93% of the spectrum delocalized, 
and the upper 7% localized, and agrees with the 
result of GNR [5,10]. 

These results make an interesting comparison to 
our earlier work on amorphous Si [11]. One major 
difference is that the density of vibrational states 
in am-Si retains much of the structure seen in the 
crystal, with acoustic and optical types of vibra- 
tions separated by a region of low density of states. 
To the contrary, in am-LJ, the density of states 
loses differentiation and is completely featureless. 
Apart from this, other properties of am-LJ vibra- 
tions are reminiscent of those in am-Si. We see sim- 
ilar behavior in the collapse of the bond-stretching 
character onto a single frequency-dependent curve. 
The anomalous states at very low frequency, and 
the strongly localized states at very high frequency, 
have similar character. In both glasses, fewer than 
10% of the modes are localized, always at the up- 
per end of the spectrum. 
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